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ABSTRACT 

We discuss the relationship of the forward matrix element of the operator '^'0, re- 
lated to the so-called sigma term, to the quark number. We show that in the naive 
quark model in the canonical formalism these quantities coincide in the limit of 
small average quark momenta. In the QCD parton model defined through light- front 
quantization this result is preserved at leading perturbative order but it receives ra- 
diative corrections. We analyze the theoretical and phenomenological consequences 
of this result, which provides a bridge between a current algebra quantity, the sigma 
term, and a deep-inelastic quantity, the parton number. 
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1. INTRODUCTION 



The fermion mass operator '^'0 plays a peculiar role in QCD, due to the fact that 
it is the only term in the Lagrangian of the theory which breaks both chiral and scale 
symmetry at the classical level. Therefore the matrix elements of this operator enter 
several low-energy relations which may be derived from current algebra using chiral 
symmetry. For instance, the vacuum-expectation value of this operator controls 
the pion mass; the nucleon expectation value of the same operator is related to the 
pion-nucleon scattering amplitudes at threshold, and so forth. On the other hand, 
the same matrix elements are related to the breaking of flavor symmetry due to 
scale noninvariance, i.e., to quark model formulas for mass splittings. Whereas the 
former set of results points towards the meaning of 'i/'V' as a probe of the low-energy 
properties of hadrons, the latter suggests that its matrix elements could be related 
to the total quark content of hadrons, a quantity which has an interpretation in 
terms of partons as seen in high-energy experiments. 

Indeed, it has been suggested that the (forward) matrix element {'tpi4>i)h 
of 'i/'i'i/'i ill a hadron h may be assumed to be proportional to the total number of 
quarks plus antiquarks of flavor i in that hadron (at least in the limit of small quark 
mass). This matrix element, summed over light flavors, is experimentally accessible 
due to its proportionality to the so-called nucleon sigma term, namely the matrix 
element of the mass term proper, which can be calculated and measured in a variety 
of ways: 

a = {h{p)\m^i)i%l}i\h{p)) =m^{il)i%l}i)h, (1.1) 

i i 

where m is the average value of the current light quark masses and the sum runs 
over light quark flavors. 

The validity of this assumption is established by its phenomenological success: 
for instance, it provides a determination of the current quark masses |jl[], and it 
turns out that the values determined in this guise are in good agreement (to about 
20%) with those found through utterly different techniques [0. This assumption 
has been subsequently used in different contexts [Q, with a varying degree of 
phenomenological success. 

Because the total number of quarks plus antiquarks in a hadron is also given by 
the first moment of the nucleon structure function F2{x) [|5[ which can be measured 
in deep-inelastic scattering experiments, if this phenomenological assumption were 
correct it could provide a very important bridge between high- and low-energy 
properties of hadrons. For example, measurements of the isotriplet sigma term from 
low-energy scattering experiments would thus provide independent information 
on the pattern of SU(2) violation in the first moment of quark structure functions 
which has been recently measured in deep-inelastic scattering experiments with 
rather surprising results. 

However, whereas the identification of {i(^i(^)h with the quark plus antiquark 
number is suggested by the naive observation that the operator t/jtl; has charge 
conjugation properties opposite to those of the charge operator tp^ip, it is clear that 
this cannot be an exact identification, because the number operator is obviously 
diagonal in quark number, whereas ipip has in general nonvanishing nondiagonal 
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matrix elements, i.e., it connects Fock states whose quark numbers differ by two 
units 0. 

Also, the link of the sigma term to a moment of a structure function cannot 
be established through the standard procedure of taking moments of the light-cone 
expansion of electromagnetic currents, given that t/jil^ has twist 3, and therefore it 
does not appear in the operator expansion at leading twist. Indeed, a more detailed 
analysis, based on the rigorous definition of parton distributions in terms of light- 
cone field operators [§, shows that the nucleon matrix element of t/^tl^ is generally 
related to a quark-quark-gluon correlation function, rather than to a bilinear quark- 
quark operator only, as the number operator is. 

It seems thus that, unless the successfulness of the identification of {ipij)h with 
the quark number [Q is a mere accident, it must be a consequence of dynamics, 
rather than symmetry. It is the purpose of this paper to investigate the relevant 
dynamical issues, in order to establish whether, and to which extent, this identifi- 
cation is correct. 

To this end, we shall study the partonic interpretation of the operator 'i/'?/', 
both in the naive parton model in the canonical formalism, and in the QCD parton 
model in the light-cone formalism. We shall see that in the naive parton model 
even though the operator '^'0 and the quark plus antiquark number differ by non- 
diagonal terms, their forward matrix elements coincide in the limit of small quark 
momenta, whereas they are proportional through a scale-dependent coefficient for 
generic quark momenta. We shall discuss to which extent this coefficient can be ap- 
proximately neglected at the nucleon scale for the experimentally observed parton 
momentum distributions. 

We shall then provide a rigorous expression of the operator •i/'?/^ in terms of the 
quark distributions defined in terms of field operators on the light cone. We shall 
show that the naive quark-model results are corrected by terms which measure 
quark-gluon correlations. These terms are superficially enhanced by a power of 
energy over mass, and are therefore usually thought to be the leading (or only) 
contribution to tpil} in the limit of small current quark mass. However, we shall see 
that the leading term in a small mass expansion vanishes, and therefore at leading 
perturbative order the naive parton model result is reproduced, whereas it receives 
radiative corrections which are higher order in the QCD coupling, while being of 
the same order in the quark mass over energy as the number operator terms. 

In sect. II we shall compute the operator ipip in terms of canonical quark creation 
and annihilation operators, we shall see that the forward matrix elements of its 
non-diagonal part in quark number vanishes, whereas the surviving diagonal part 
coincides with the quark plus antiquark number operator up to a scale-dependent 
constant, and we shall estimate the magnitude and mass-dependence of this constant 
using well-established parametrizations of the quark momentum distribution. In 
sect. Ill we shall relate the forward matrix elements of the operator ipil) to the 
parton distribution functions as defined using the QCD parton model in the light- 
cone formalism, and we shall see that at leading perturbative order the naive quark 
model result is reproduced. We shall conclude with a discussion of our results in 
sect.IV. 
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2. THE OPERATOR i/^?/^ IN THE CANONICAL FORMALISM 



Use of the matrix elements {il'il^)h of the operator t/ji/; as a measure of the 
number of quarks plus antiquarks in the hadronic state h is suggested by the naive 
observation that the operator ^fjifj has charge- conjugation properties opposite to 
those of the charge density i/j^i^, which coincides with the difference of the number 
operators for quarks and antiquarks. In order to derive quantitative consequences 
of this simple remark, we compute the operator -i/'-i/' in terms of quark and antiquark 
creation and annihilation operators. 

To this purpose, we introduce the standard expansion of the fermion field 



(27r)3/2 



1/2 



hr{k)ur{k)e-'^'' + dl{k)vr{ky 



[2.1) 



where, in general, the operators 6^^, d) create quarks and antiquarks, respectively, 
in the interaction representation, and ij} reduces to a free field operator only at 
asymptotic (initial and final) times 0. Since we are interested eventually in the 
matrix elements of ipil} at zero momentum transfer, we compute its integral over all 
space, with the result 



+ 



s 

+ E/^'^| [M.,.(^)+Mt,(fc) 

iV(+)(fc) = h\{k)h,{k)- NtHk) = dl{k)ds{k) 
Mr,s{k) = h\(k)dl{-k) \us{k)vr{-k)e^' 



iEt 



(2.2) 

(2.3) 
(2.4) 



where the operator has been canonically normal-ordered. 

Eq. (|2.2| ) shows that at zero momentum transfer decomposes in the sum 
of four distinct contributions. Two of these are diagonal in quark number, and 
contain the canonical number operators for quarks and antiquarks of momentum 



k and spin s, Ns^'{k) and Ns \k), Eq.( |2.3|) , whose matrix elements are simply 
the number densities of quarks or antiquarks with the given spin and momentum. 
The remaining contributions, instead, create or annihilate qq pairs. For example, 
Eq.( p.4|) creates a quark and an antiquark with momenta k and —k and spin r and s 
respectively. All of these operators are in general time-dependent; the diagonal ones 
( p.3| ) implicitly through the time-dependence of the interaction-picture operators b 
and d, and the non-diagonal ones ( p.4|) both implicitly, and explicitly due to the 
energy exponential. 

Because the operator ifjifj depends generally on time, its physically measured 
matrix elements should be defined as time averages: 



lim — 

T — >oo 7" 



72 



72 



dt {h\ / d^x -.tpij: \h), 



(2.5) 
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where r is a proper time interval, so that the matrix element is a Lorentz scalar 
as the operator t/jt/; is. This definition coincides, up to the cofficient m, with the 
standard definition of the a term ( p. . 1|) , as it is clear by noting that they are both 
Lorentz scalars, and they coincide in the rest frame of the state \h). 

The matrix elements of the nondiagonal operators ( |2.4| ), instead, seem to spoil 
the claimed identification of 'i/^'i/' with the quark number. We shall now show, how- 
ever, that upon averaging according to Eq.( p.5|) these matrix elements vanish, with 
mild assumptions on the state \h) in which the average is taken. Introducing the 
usual u and v spinors, namely 



"•« = \/^(^)^'' (2.6) 



Vs{k) = Z7 ul{k), 

where Xs is the two-component spinor of spin s, the non-diagonal term of Eq. (|2.2| ) 
can be written as 

^ = E/ f Mr,s (^) = - E / ^e'^^'Hmli-k) [xla ■ kx-r] , (2.7) 

r,s r,s 

plus its hermitian conjugate which we shall drop henceforth. Performing the 
spin sums explicitly obtains 

-^e^'^'[bl{k)dl{-k) + bl{k)dU-k)]k, 

_ jd^ e^^Et^bl(^k)dl{-k) - bl{k)dl{-k)]k^ (2.8) 

/d^ k 
— e^''''[bl{k)dl{-k) + bl(k)dl{-k)]ky. 

The various operators which appear in Eq. (|2.8| ) can be identified with those 
which create a qq pair with definite spin: 



(2.9) 



These are coupled to angular distributions corresponding to L = 1 states 
L = l = k, = k^Y^{eA) 

L = l L^ = -l -k^ + iky = -k.f^Y-\e,(j)) (2.10) 



s = 


1 


s, = o 


b^ d^ + 


s = 


1 


s, = i 


b\d\ 


s = 


1 


s. = -i 


bldl 



L = l L^ = l k^ + iky = -kJ^Y^{e, 
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so that the operators in B create the J = L + S = state 



Y k [5'i,iLi _i - Si^oLi^Q + 5*1 _iLi,i] = VStt/c Jo,o, (2.11) 

where by 5'i,i we mean a state with S = 1, Sz = 1, and so forth, and 9, (p and r are 
the spherical components of the q — q relative three-momentum k. 

Consider now the time average of the matrix element of the operator B Eq. ( |2.7| ) 
in a state \h). If this were just a free- field state, then the matrix element would 
vanish trivially because of energy conservation, since upon integration over t the 
energy exponential yields a Dirac delta which in free field theory would force the 
energy of the pair to vanish. However, with the application to hadronic states in 
mind, we may assume the state \h) to be a Fock space superposition of many- 
particle bound states. In particular, we assume the state \h) to contain couples 
of Fock states which carry the same charge but differ by two units in quark plus 
antiquark number, and such that the extra q-q pair forms a bound state described 
by a wave function '(/'(/c). The (forward) matrix element of B in this state is then 
(according to Eqs. (|2.9|) -( |2.11| )) proportional to the J = state wave function, with 
a coefficient C2 depending on the overlap of the two Fock states once the pair is 
removed. This, in general, may be nonzero even after imposing energy conservation, 
because of the presence of a nonzero binding energy. 

Assuming the state \h) to be unpolarized, the matrix element of B in the 
hadron's rest frame is thus given by 

{h\B\h) =C!,J ^e^'^'^kRik) (2.12) 

where R{k) is the radial part of the q-q wave function, we have absorbed dimen- 
sionless factors in the coefficient C2, and E' = E — Eq = \/k'^ -\- m? — Eq, where 
2Eq is the binding energy of the bound state. After time averaging according to 
Eq.( ^75|) we get 

{B)h = C'^ lira - r dt I — e^'^''kR{k) = (2.13) 
T^oo r J E 

1 k^ 

= 47rC2 lim - / dk —— R{k) sin E't. (2.14) 

T^oo T Jq EE 

It is easy to see that assuming a simple Gaussian radial wave function 

R{k) ~ e ^, (2.15) 



the momentum integral in Eq.(2.13) 



»oo 7.3 

lit) =4n I dk^ R{k) e^'^ ' (2.16) 
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satisfies lim,-^oo 7 J-rdt I{t) = 0, since J^^dt I{t) is finite. Tliis proves tliat, with 

this choice of radial wave function R{k), the limit ( |2.13|) is zero. The argument can 
be generalized to anylil radial wave function such that the average value of k"^ 

/■OO /"OO 

{k^)= dkk^R^{k)= dE E{E'^ - m^)^ R^{E) (2.17) 

Jo Jm 

is finite. Indeed, demanding a finite value of (/c^), Eq. ( |2.17| ), implies 



R{E) ~ E-i-^ 
R{E) ~ {E-m) 



(2.18) 



SO that the function ^-^P^R[E) is integrable both when E m and E 00. It 
follows that 

1 /"°° k^ 
{B)h = 4tvC2 lim - dk— — R{k) sin E'r = 
T^oo r Jo EE 

= 47rC2 lim - / dE R{E) sm{E - Eo)t = 



(2.19) 



where the last step is a consequence of the Riemann-Lebesgue lemma. 

It follows that the time average of the matrix element of the non-diagonal 
contribution Eq. (|2.4|) to ipip vanishes under the assumption that the nondiagonal 
operator creates (or annihilates) a pair describing a bound state with finite average 
momentum. If we consider now the forward matrix element of ?/''!/' in a hadronic 
state in the center-of-mass frame, we may view in a naive parton model approach 
the states created and annihilated by tj^ip as quark-parton constituents, and this 
assumption translates into a very natural assumption of boundedness of the quark 
momentum distribution, which is satisfied by phenomenological quark-parton wave 
function considered in the literature [l^]. We conclude that the hadronic matrix 
elements of 'i/'-i/', Eq.(|2.5|), reduce to those of its diagonal portion: 



/^^)^ = y lim - [ ^ dt(h\ [d'^k- \Ni+\k) + Ni-\k)\ \h). (2.20) 

A^r^^ T J_^/2 J E I J 



r/2 

I 

In the naive quark model Eq.( p.2Ci|) simplifies further to 



W). = ^ l^[q{x)+q{x)]. (2.21) 



^ We assume for simplicity R{k) to be real; this does not entail loss of generality 
because, according to Eq.s ( |2.2D , ( |2.7| ), the matrix element of B contributes to {ipip) only 
in the combination B + BK 
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where q{x) is the quark distribution and m and M denote the quark and 



hadron mass, respectively. In Eq. (|2.21|) manifest Lorentz invariance, which was not 
apparent in the intermediate steps of our computation, has been recovered. This 
Eq. provides the desired expression of the forward matrix element terms 
of quark distributions in the naive quark model. 

The quark plus antiquark number is given instead by the first moment of q{x) 

n = jdx [q{x) + q{x)] , (2.22) 

which diverges in the singlet case, but is finite if one takes nonsinglet matrix elements 



T0[] . Hence, the diagonal portion of '^'0 differs from the sum of canonical quark 
and antiquark number operators due to the presence of a coefficient of ^, which is 
necessarily present because ip^ip is a Lorentz scalar (as opposed to the charge density, 
which is a component of a four- vector). This coefficient makes the identification of 
the operator ipip with the quark number rather puzzling in two different respects. 
On the one hand, for a given quark fiavor, this coefficient depends on the energy 
distribution of quarks, and differs significantly from unity if the distribution is 
hard, i.e., if it contains a significant fraction of large-momentum quarks. On the 
other hand, the coefficient depends on the quark mass and is thus highly fiavor- 
asymmetric. 

Due to the first effect, the matrix element of 'i/'?/' in a hadronic state is at best 
proportional to the quark content, rather than equal to it: 

= C^n,, (2.23) 

where the coefficient of proportionality is roughly given by = (^) the average 
value of the coefficient for the given fiavor. Due to the second one, this coefficient 
depends on the mass. This means that if one defines, as in Ref. ||l|], the current 
masses as the coefficient of proportionality between the sigma term Eq.( |l.lj ) and 
the quark number, then the values of the masses thus defined are related to the true 
values (defined, for example, as chiral symmetry breaking parameters 0) through a 
coefficient that differs from unity and depends on the mass itself. This is in striking 
contradiction with the phenomenological observation that the values of the masses 
determined in this guise agree to about 20 % with those determined using chiral 
symmetry, thereby suggesting that this coefficient is compatible with unity, and 
moreover mass-independent. 

Because the quark momentum distribution is at present uncalculable, this 
should be taken as a dynamical accident. We may however understand the ori- 
gin of this effect by considering the structure of quark distributions. These can be 
phenomeno logically expressed as the sum of a valence part, whose first moment is 
finite and equal to the naive quark content of the hadron, and a sea part, whose 



first moment diverges. Valence-quark distribution behave as q{x) ~ |]TT|], con- 



sistently with the idea that valence distributions stem from mesonic Regge pole 
exchange. This implies that if we neglect sea contributions (having the application 



to nonsinglet matrix elements in mind) then the integral on the r.h.s. of Eq.( 2.21 
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diverges at small x as ~ X-^^.^. Now, in the naive parton model the range of values 
that X may take is bounded dynamically (by the requirement that partons carry a 
finite fraction of the parent hadron's momentum), leading to the value Xmin = (j^) 
p2} j. It follows that the integral on the r.h.s. of Eq.( |2.2lD , which is dominated by 
the small-x region, behaves as 

m f dx , . m -I /m\0 , 

M j -«'^l~aM^"'".-o (m) • (2-2*' 

i.e., it does not depend on the current mass m. Thus, if we consider nonsinglet 
matrix elements of 'i/'?/', and we assume that the sea contributions cancel, we see 
that, contrary to naive expectation, the coefficient of ^ in Eq.( p.20| ) does not lead 
to any dependence on the quark mass. 

Actually, the assumption that the nonsinglet sea quark distributions should 
vanish has been recently falsified experimentally 0,0. However, Regge theory 
leads to expect [|TB[ that the nonsinglet part of the sea distributions behaves at 
small X in the same way as the valence distributions; this is borne out by fits to the 
available data |13]. It follows that the result ( p.24|) holds even if this effect is taken 



into account. 

The divergence of the valence distribution at small x, besides explaining [due 
to Eq.( |2.24| )] the apparent mass-independence of the coefficient Ci [Eq. ( f2.23| )] also 



shows why this coefficient is close to unity. Indeed, this divergence implies that in 
the matrix elements of ^/^V the integration over quark momenta is dominated by 

2 

momenta of order k ~ Mxmin, i.e., k ^ ^ << m (for the three light fiavors), so 
that f ^ 1. 

In sum, the valence quark distributions appear to be soft enough that the coef- 
ficient Ci in Eq.( p.23|) may be approximated with unity, and their small-x behavior 
is such that this coefficient is mass- independent. This, together with the vanishing 
of the matrix elements of the non-diagonal portion of 'i/''^, provides a phenomeno- 
logical justification for its identification with the quark number operator in the 
naive parton model, at least for the three lightest fiavors (notice that in Ref. [0 this 
identification was suggested in the m — limit). Even though this result is intu- 
itively appealing, it does not illuminate the relation of to the parton operators 
which can be constructed in QCD, and related rigorously to the matrix elements 
which appear in the leading order of the light cone expansion which is relevant for 
deep-inelastic scattering. In order to do this, we must go from the canonical, naive 
parton model to the QCD parton model on the light cone. We shall do this in the 
next section. 
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3. AND THE LIGHT-CONE QUARK OPERATORS 



In the previous section we derived a relation between the operator ipip and 
the quark number within the naive parton model picture of a nucleon made of 
quasi-free constituents which share its momentum, which were identified with the 
states created by the interaction-picture coefficients of the Fourier modes of the 
quark field. However, in deep-inelastic scattering the nucleon's constituents are 
only probed through the measurement of structure functions. The moments of 
structure functions are then linked to the nucleon matrix elements of the twist- 
two operators which appear in the leading order of the light cone expansion of two 
electromagnetic currents. These are related to physical observables (such as the 
quark and gluon number or momentum) that can be expressed in terms of quark 
and gluon distributions, which are the basic quantities in the QCD parton model 
1^. Because the operator ip'ip is twist-three, it does not relate directly to one such 
observable, and it may be expressed in terms of properties of the quark and gluon 
distributions only by means of its interpretation in terms of quark and antiquark 
field operators. 

In order to establish a rigorous relation between the matrix elements of il^ij), 
and the quark distributions whose moments are measured in deep-inelastic scatter- 
ing we need therefore a formalism which allows to express the quark field operators 
in terms of the quark distributions themselves. Such a formalism is provided by 
the light-cone approach P], where field operators are quantized on a light front. 
This allows a separation of the physical quark degrees of freedom, by means of a 
decomposition of the fermion field operators into unconstrained components, which 
describe the dynamical degrees of freedom (called "good" components), and con- 
strained components, expressed in terms of the former and of the gluon fields by a 
constraint equation ("bad" components). We shall follow Ref.Q, to which we refer 
for a detailed treatment; notations and conventions are summarized in Appendix 
A. 

In the light cone approach quark-parton creation and annihilation operators 
are defined by the Fourier decomposition of the good components of the Dirac field 
at y+ = 0: 

1 dk^ C ^-tr-^ r ^ ^1 

V+(0,^,2/-) = ^ I ^ jdkY,[Uk^.mnse-'^y + dl{k+ , k,0)v, e^'y , 

(3.1) 

where ky = k~^y~ — k- y. The quark-parton distributions of a state |P) are then 
simply given by 

■)(P|^;(0,y,y-)^+(0,0,0)|P) 

^Hr(P|^+(0,y,2/-)^;(0,0,0)|P) 

(3.2) 



'.(x,^l = ^^/c^y-/c^ye- 
'd-rk) = ^^Jdy-ldye-^ 
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g(x) = -^ / dy- e-^^^'^y- {P\i;l{0,6,y-)M0,6,0)\P) 
q{x) = ^ J dy-e-'^^^y~tv{P\^+{0,6,y-)^l{0,6,0)\P), 



where |P) are momentum eigenstates in coordinate space, normalized as {P\P') = 
{27T)^2P~^dss'S{P~^ — P~^')d^'^\P — P'). Upon integration over the transverse mo- 
menta k Eq.s yield 



(3.3) 



In terms of the creation and annihilation operators defined in Eg .( 13 .11 ) the quark- 
parton distributions read 

= aAsTT I dky2{h\bl{xP+,k,0)bs{xP+,k,0)\h) 
dkJ2{hWi^Hk)\h) = jdkVq{x,k) 

1 1 r ^ ^^-^^ 

^(^) = dky2{h\dl{xP+,k,0)ds{xP+,k,0)\h) 
(27r)"^ 2x J ^ 

dkJ2{hWt\k)\h) = [dkVq{x,k), 



1 1 

f27r)3 2^ 



1 1 



{2ny 2x 

where \h) are momentum eigenstates in momentum space, normalized to {h\h) = 1. 
This equation expresses the quark and antiquark distributions q{x) and q{x) as the 
forward matrix elements of the number operators A^i^'' (k) and Ng (k) which count 
the number of quarks and respectively antiquarks with momentum k and spin s in 
the given state at fixed "time" y~^ = 0. 

The moments of quark distributions defined by Eq.s(|3.3|), (|3.4p satisfy the usual 
relations to matrix elements of twist-two operators . Higher- twist operators, such 
as '^'0, are not given by matrix elements of the good components of the quark field, 
but rather, they contain both good and bad components; each unit of twist carries 
an extra bad component. Since bad components are functions of the good ones and 
of the gluon fields through the constraint equation that they satisfy (see Appendix 
A) , higher twist operators describe generally multiparton correlations [jl^ . 

In the particular case of the operator we are interested in we have 

= i= ij\P.^+P+ + P+^P-)^ = (3.5) 
= -^^\l~tp- +h.c., 
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which in terms of good component only reads 



v2 



(3.6) 



where the suffix ± labels the transverse components (1,2). This is the expression 
which we ought to analyze in terms of quark (and gluon) creation and annihilation 
operators in order to understand the rigorous relation of '^'0 to the quark number 
(if any). 

The forward matrix element of 'i/'i/' can be defined as the first moment of the 
light-cone Fourier transform of the combination of quark bilinears appearing in 
Eq. ( p75|) , analogously to what was done in Eq. (|3.3| ) for quark and antiquark distri- 
butionsB 



(PIV^O, 0, y-)7-^-(0, 0, 0) + h.c.\P) + 



+ tr(P|i/;_(0, 0, 2/")i/'i(0, 0, 0)7- + h.c.\P) 



^jdx dx {h\Aic + Bic\h) , 



(3.7) 



where in the last step we used Eq. 



and 



Ale = -im :'i/'+(x)(a+) ^ijj+ix): 



(3.8) 



and 



(3.9) 



Notice that whereas Aic is a bilinear quark operator, Bic correlates a quark bi- 
linear and the gluon field operator j^^. The expression Eq.(p77|) of {^'^)h can be 
constructed to be manifestly a Lorentz scalar by averaging with respect to y"*", as 
discussed in the previous section (Eq. (|2.5| )). 



^ The normalization is fixed by defining the generic bilinear '{[jVip and requiring that 
when r = \/27^ then the first moment be equal to the total quark charge. The factor 
of \/2 is due to the requirement that in the nucleon's rest frame the second moment of 
the quark minus antiquark distribution be equal to the hadron's mass M = y/2P+. The 
matrix element of tptp is then obtained by setting F = 1 . 
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Let us now compute the contributions from terms (^]^) and (|3.9| ) separately. 
Using the Fourier decomposition of the field operator Eq.( |3.1D we get 



dx-dx{h\Mh)=^^jdx-dxj^ St/ ^rjd^jd^' 

s,s' ^ ' 

+ h\{k)d\,{k')e'^''+'''^^ ul{k)vAk')-^+ (3-10) 

iik^) 

+ ds{k)hAk')e-'^^+'''> vl{k)uAk')TAr 



{-ik'+) 

d\,{k')d,{k)e-^^^-^')- vl{k)vAk')jr^Ah). 



[ik' 

where bl{k) stands for 6^(/c+,/c,0) and so forth. The terms proportional to the 
operators b^d^ or db which create or annihilate a quark-antiquark pair in the above 
equation vanish, since the dx~ integration provides a d{k~^ + k'~^) which yields zero 
upon integration over dk'~^ . The light-front quantization thus eliminates automati- 
cally the spurious contributions to i/jiIj which vanish upon time-averaging, discussed 
in the previous section, and corresponding to pair creation, in that the physical 
quark degrees of freedom are selected by the projection of the good components of 
the Dirac field. 

Exploiting the explicit expression of the Dirac spinors given in the Appendix 
we remain with 



dx dx {h\Aic\h) — 

°° dk+ r dk 



(3.11) 



= vi^Jo + 4^(^)1/.). 

Setting = xP~^ we get finally 

/TTl f dx 

dx- dx {h\Ai,\h) = J — [qix)+qix)], (3.12) 



where we have used the definition ( |3.4| ) of the quark densities. We may now exploit 
again the fact that {^ip)h is a scalar, and evaluate the r.h.s. of Eq.( |3.1'^ ) in the 
hadron's rest frame, where = The r.h.s. of Eq.( p.l2D is then recognized to 



coincide with the expression of the operator 'i/'i/' Eq.s( |2.2D ,( |2.21| ), which we derived 
in the previous section in the naive parton model. 

It follows that the purely fermionic contribution to the operator Tpip reproduces 
the naive parton model result. The full matrix element of ipip, however, reduces to 
this only in the limit of vanishing quark-gluon coupling (7 — 0, since the Lorentz 
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transformation required to transform to the light-front the interaction picture cre- 
ation and annihilation operators has brought in a dependence on the interaction, in 
the form of the terms Bic, which depend on the gluonic degrees of freedom. 

Superficially these terms provide actually the dominant contribution to the 
matrix elements of ipip in the small mass limit, in that they do not carry an explicit 
mass prefactor; therefore the operator '^'0 is sometimes identified with a quark- 
gluon correlation, and its portion related to the quark number, Eq. (|3.12| ), is entirely 



neglected This would of course spoil the identification of ipip with the quark 
number altogether. Let us thus turn to the contributions of terms Bic Eq.( |3.9| ), in 
order to see whether this is really the case. To this purpose, we need the Fourier 
expansion of the gluon field operator 



1 



dp 



+ 



2p^ 



dp ax{p)e^^e-'P^ + a\{p)t^^{p)e^P- 



(3.13) 



where the operator a\{o}-^ annihilates (creates) a gluon with four-momentum p^ 
and polarization vectors (see the Appendix). 

We compute Bic, putting the Fourier expansions of the quark and gluon fields, 
Eq. (pT|) and (|3.13| ), in its definition Eq. (p^) , with the result 



j dx dx :'i/'+(a^) '^■:4-a^ 



dkdk' 



dk+ r dk+' 



-I 1 



+ 



2k+ 



°° dk+ dk+' 



1 



2k+ Jo 2k+^''^''^'^'^'''^''^^^~^'^'^^'^^^''''^^^Jl^-k+')^ 



+ h.c. 



(3.14) 



The matrix elements of this operator may now be determined explicitly using 



u 



u 



:{k) /I uAk') = V2k+k+'ds,-s'[2s5x,i - i5x,2] 



Xk) Vs^{k') = V2k+k+'ds,s'['^s6x,i - t6x,2] 



(3.15) 
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which leads to 



dx dx {h\Bic\h) 



9 



2v^(2 



dkdk'V2k+k+' 



X {hi 



dk+ r''^ dk+ 



bl{k)b-{k')a+{k - k') - bl{k)b+{k')a-{k - k') 



,0 2k+ Jo 2k+' 



+ 



dk+ 



dk+' 



{k+ - k+'Y 

b\{k)d\{k')a+{k + k') - bl{k)dl{k')a-{k + k') 



2k+ Jo 2k+' 



{k+ + k+')^ 

dk+ f^^ dk+' [dl{k)d-{k')a+{k - k') - dl{k)d+{k')a-{k - k')] 



2k+ Jo 2k+' 



{k+ - k+') 



i\2 



+ h.c. )\h), 



(3.16) 

where o}±{p) = -^[a\{p)±a\{p)] creates (and a± annihilates) gluons with ± helicity. 

Eg . (13.161) provides the general form of the quark-gluon correlation measured 
by the matrix elements of '^'0, and shows that Bi(. connects Fock components of the 
hadron wave function which differ by one unit in particle number. More specifically, 
there are three classes of contributions to this matrix element: (a) correlations of 
a quark state to a quark-gluon state (or vice- versa); (b) correlations of a quark- 
antiquark state to a gluon state (or vice- versa); (c) correlations of an antiquark 
state to an antiquark-gluon state. Due to the helicity structure of the operator 
i/'i/j, in all of these the fermion helicity is fiipped at the vertex corresponding to the 
insertion of '^'0, i.e., to the effective interaction Eq. (|3.9|) , while the gluon carries 
away one unit of helicity, thereby ensuring angular momentum conservation. 

The matrix element, Eq.( pj.l6D , is therefore given by the various overlaps which 
measure the amplitude for finding the pertinent quark-gluon states in the hadron's 
Fock wave function: 



{h\Bi,\h) = Ba + Bb + B, 

Ba = {h\qg){qg\BiM){q\h) + {h\q){q\Bi,\qg){qg\h) 
Bt = {h\g){g\BiMq){qq\h) + {h\qq) {qq\Bic\g) {g\h) 
Be = {h\qg){qg\BiM){q\h) + {h\q){q\BiMg){qg\h), 



(3.17) 



where Ba, Bb, and Be correspond respectively to the three classes of contributions 
mentioned above. 

All the matrix elements in Eq.( |3.17] ) can be computed in terms of the quark 
distributions of the given hadronic state \h), if we make the assumption that the 
gluons and antiquarks are generated by QCD radiation processes from the quarks. 
This assumption is eminently reasonable in the light-cone formalism, where [Q] the 
dynamics can be thought of as that which takes place in the reference frame where 
the state \h) moves at a speed tending to the speed of light (the formalism is 
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constructed to be invariant upon the required boost) and therefore the constituents 
of \h) are quasi-free. 

Let us consider for definiteness the matrix element S^, Eq.( |3.17|) . We get 

Ba = {h\q){q\C,nt\q9){q9\BiM){q\h) + {h\q){q\Bl,\qg){qg\C^nt\q){q\h), (3.18) 

where 

C^nt =9 4 V^: (3.19) 

is the QCD interaction Lagrangian. Similar expressions may be derived for the 
other contributions Bi, and B^. to the matrix element of S/c- 

We may work out explicitly the various matrix elements in Eq.( 3.17] ) by spelling 
out the interaction lagrangian (|3.19[ ) in terms of quark and gluon creation and anni- 
hilation operators, analogously to what was done for the operator Bic in Eq.( |3.14] ). 
We get 



= / dkdk' / dy+ J] 

^ > -J J s,s',X 



dk+ r'^^ dk'+ X p*^'' ~p ^y'^ 

y^blik)bAk')ax{k - k')usik) j. iX)uAk') 



2k+ Jo 2k'+ «^ ^ ^ ^ ^ y .V . r V . . V . k+ - k'+ 

dk+ r°° dk'+ X 4 pi{k-+k'--p-)y+ 
+ / ^/ ^hl{k)dl{k')a,^{k + k')u,{k)n^)vAk'y 



2k+jQ 2k'+ ' ^'^ ' . .V .rv . .V . k+ + k'+ 



dk+ f^'' dk'+ . 



-dl,{k)ds{k')ax{k - k')vs{k') ^ {X)v,'{k)- 



J(k -k' -p )y 



+ 



2k+jQ 2k'+ ^ ^ / .V /rv / / k+-k'+ 

+ h.c. 

(3.20) 

The terms in this expansion provide the transitions corresponding to the various 
QCD radiation processes: for instance, the first term corresponds to gluon absorp- 
tion by a quark, the second to a pair creation by a gluon, and so forth. 

The computation of Ba reduces thus to the computation of all the Feyn- 
man diagrams which contribute to the transitions {q\'Cint\qg){qg\Bic\q) and 
{q\Bic\qg) {qg\Cint\q) ■ These are displayed in Tab. 1, along with the result of the 
computation. Inspection of these result indicates that the contributions to Ba fall 
in two subclasses, namely, those where the helicity of the initial and final quark is 
the same, and those where the helicity is changed. The former are suppressed by 
a power of mass (that is, ■^) with respect to the latter, due to the fact that they 
require a helicity fiip in the quark line at the perturbative gluon emission vertex. 
It follows that the leading contributions in the limit of small quark mass are those 
which correspond to interference between different helicity components of the given 
state \h), i.e., those which are proportional to the density of transversely polarized 
quarks in the state. 
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We may now determine Ba explicitly by using in its expression Eq. (|3.18|) the 
results listed in Tab. 1. Notice that care must be taken in the first (second) term 
of Eq. (|3.18| ) to ensure that the QCD interaction takes place, respectively, at a 
later or earlier "time" compared with the x"*" = "time" of the interaction ( |3.14|) ; 
this is accomplished by introducing suitable Heaviside 9 functions. A tedious but 
straightforward computation leads to 



V2mg^ f^dk+ [>• dk+' , J J, 



_ _ , (3.21) 

-i(k —k' —p 



X 



/ -i(fe -k -p ^ 



where Vq{x,k) has been defined in Eqs. (|3.2|) -(|3^). Notice that Vq{x,k) has the 
dimensions of an inverse mass squared, due to the normalization of the creation 
and annihilation operators (see the Appendix). Remarkably, it turns out that all 
contributions where the quark's helicity is changed, which would be leading in the 
small mass limit, cancel against each other. 

A similar computation leads to the same result for the terms B^ with Vq{x, k) 
replaced by Vq{x,k), the probability of finding an antiquark in the given state. 
Finally, the terms Bi, are equal to 

(3.22) 

where Vg{x,k) is the gluon distribution, defined in analogy to Eqs.(p72| ) - (|3.4|) . It 
follows that all the leading contributions in the small mass limit to Bic vanish; the 
only surviving contributions are those suppressed by a mass prefactor, i.e., those 
which are of the same order in mass as the operator Ad, Eq. (|3.8| ), proportional to 
the quark plus antiquark number. 

Using the integral representation of the Heaviside function 



e{y+) = -2ik'^e V '"'^ J / ^ ^ (3.23) 

J_oo 27r k' — m? — ie 

in Eqs. (|3.21| ) and ( |3.22| ) allows to perform the integration, with the result 



{h\Bi,\h) = Ba + Bh + B^ 
m f dx 



V2P+ J X 



— dk 



Vq{x, k) + Vq{x, k) ) gi{x, k) + Vg{x, k)g2{x, k) 

(3.24) 
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where the functions gi are exphcitly given by 

,2 /-xP"" 



gJx,k) = ^^xP+ dk'+dk' 



1 



X 



k'+{k+ - k'+){k^ + m2) - k+k'+{k - k'f - k+{k+ - k'+){k'^ + m?) 

g2(x,k) = —^-^XP+ dk'+dk' 
TT^ 47r Jo 

1 



2k+k'+k ■ k' - k'+'^k^ - k+'^k'^ - Ak+'^m? 

(3.25) 

2 

and after renormahzation of the loop integration over /c', is replaced by the 
running coupling ctg. 

Hence, explicit computation reveals that the leading contribution to the portion 
of which is proportional to a quark-gluon correlation, Bic Eq. (|3.9| ), vanishes, 
leaving a term which is of the same order in the quark mass as the terms Aic, related 
to the quark plus antiquark number according to Eq.( |3.12]) . This surviving term 

is proportional to the quark, antiquark and gluon distribution functions V{x,k), 
defined in Eq.s (|3.2|) -(|3.4|); it is, however, higher order in the strong coupling ag- 
Putting Eq.s and (|]2|) together we get thus 




Vqix,k)+Vqix,k)] l + gi{x,k) 
+ Vg{x,k)g2{x,k) y 



(3.26) 



It follows that the naive parton model result, Eq.s( p.20D ,( p.21| ), is reproduced up to 
higher order perturbative corrections (given by the functions gi and g2 in Eq. (|3.26]) ). 
These are of the same order in ^ as the tree-level result. 

Even though we cannot determine the magnitude of these corrections, since this 
would require knowledge of the k dependence of the functions V, we may determine 
the small-x behavior of the functions gi Eq. (|3.25| ), in order to verify that the small-x 
behavior of the tree-level result is not spoiled by the corrections. This behavior, 
according to Eg. ( p. 24 ) is responsible for the independence of the quark mass of the 
coefficient which relates the matrix element of t/jtl^ to the quark number. 

The terin containing g2 in Eq. (|3.26| ) does not contribute to nonsinglet matrix 
elements of ?/''!/', and the small-x behavior of gi is 



gi{x,k) ~ 

X— >0 



dk'- 



1 



In 



P+ -ai 



In^ 



27r2 47r ./ V"^ ^+ - ^2 as 

2 



{m^ + k-k')^ -{k^ +m^){k' +m^) 
w? + k-k' ± y/A 



(3.27) 



0-1,2 



]^2 j^2 
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i.e., the function gi{x) tends to a constant at small x. It follows that the small-a; 
behavior of the integrand in Eq. (|3.24| ) in the nonsinglet case is dominated by the 
small- a: divergence of q{x) and is thus the same as that of the naive result, Eq.( |2.24] ), 
which is thus not spoiled by the first-order perturbative correction. 

We have thus shown that even though in the QCD parton model the naive quark 
model result Eq.( |2.21| ) receives radiative corrections, these, at one-loop order, are of 
the same order in ^ as the tree-level result, and do not spoil its small-x behavior, 
which controls the approximate identification of {^ip)h with the quark content of 
the hadron h. Carrying this check through higher perturbative orders would require 
a higher-order determination of the various overlaps which appear in Eq.f p.lS]) . 



4. DISCUSSION 

In this paper, we have discussed the interpretation of the fermion mass operator 
'^'0 in terms of quark field operators. We have shown that in the naive parton model 
the forward matrix element of this operator reduces to a quantity which is closely 
related to the quark plus antiquark number, Eq.s. (|2.20| ), (|2.21| ), up to a constant 



which we have argued to be mass independent and roughly close to unity. If the 
quark and antiquark numbers are defined as the first moment of the quark and 
antiquark distributions constructed rigorously in terms of light-cone parton field 
operators this result remains true at leading perturbative order. 

The origin of this result can be understood by considering the different role 
that the operator ipij) plays in the short distance expansion at fixed time, or in the 
light-cone expansion at fixed x"*", of the operator-product of two electromagnetic 
currents. In the short-distance expansion operators of low dimension dominate, 
and appears at leading order; in the light-cone expansion operators of low twist 
dominate, and only appears at next-to-leading order. It is therefore to be 
expected that this operator has a simple interpretation in terms of canonical field 
operators, and a more contrived one in terms of light-cone field operators. Since 
only the latter are relevant to partons as measured in deep-inelastic scattering, this 
would seem to spoil the possibility of endowing -i/'-i/' with a partonic meaning, other 
than as a multiparton correlation. 

Surprisingly, we have seen however that under the reasonable assumption that 
the gluon parton content is generated radiatively by the quark partons, the con- 
tribution to the matrix elements of ipil) which is leading in a small mass expansion 
vanishes. The remaining contribution is proportional to the quark's current mass, 
and is therefore related to the interpretation of the matrix element of '^'0 as a chiral 
symmetry breaking parameter. This, to leading perturbative order, is the number 
of quarks plus antiquarks, up to a coefficient of ^, where E is the quark energy. 

This coefficient is necessarily scale dependent, as it is demonstrated by the fact 
that, because the combination is renormalization- group invariant, i/^'i/' evolves 
as the quark mass, i.e., it evolves at one loop in QCD 0. The quark number, 
instead, evolves only at two loops 0; therefore the coefficient which relates and 
the number at the operator level must evolve at one-loop, too. Physically, this is 
understood as reflecting the fact that the coefficient which relates to the number 
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operator, Eq. (|2.21| ), is essentially the average quark energy, which also evolves at 
one loop. 

Phenomenological results obtained identifying the matrix elements of i/ji/j with 
the quark number at the nucleon scale suggest however that this coefficient 
should be close to unity and mass-independent. Surprisingly, assuming the small- 
X behavior of structure functions which Regge behavior leads to expect [rU|, ||ll|| 



is enough to ensure mass independence, Eq. (|2.24| ). Also, the divergence of quark 



distributions at small x justifies the proximity of this coefficient to unity. Even 
though the matrix elements of t/jtl^ are in practice only measurable at the nucleon 
scale, using current algebra 0], the approximation of these matrix elements to the 
quark plus antiquark number should actually get better at larger scales, where quark 
distributions become more and more peaked at small x. 

Finally, it is interesting to notice that the (negative) moment of quark distribu- 
tions to which i/ji/j is related, Eq.s( 2.21 ),( |2.24| ), has a divergent behavior at small-x. 



according to current knowledge of quark distributions. Nevertheless, all the result 
derived in the present paper apply to the parton model with finite quark mass, 
where the range of x is bounded from below by Xmin = (^) [l^], thus ensuring the 



mass independence of the moment to which is related, according to Eq.( f2.24] ). 
This suggests that actually the chiral limit of the forward matrix element of i/jiIj 
is smooth, a conclusion which agrees with the experimental fact that the nucleon 
sigma term, related by a simple coefficient of proportionality to this matrix element, 
is numerically close to present determinations of the nonsinglet quark number 0, 
and the theoretical fact that the sigma term is well defined and smooth in the chiral 
limit 0. 

In sum, the forward matrix element of is providing us with a remarkable 
connection between low energy properties of nucleons, as measured by current al- 
gebra relations, and their high-energy properties, as manifested by their partonic 
content. Recent precision measurements of structure functions 0, especially in the 
small- a: region, have opened the possibility of using this relation as a means to gain 
insight into parton content of the nucleon 0. The forthcoming experimental re- 
sults on the small-x behavior of structure functions hold the promise of allowing, in 
reverse, to use high-energy data as a means to uncover the structure of the nucleon 
beyond perturbation heory. 

Acknowledgements: We thank E. Predazzi for discussions and D. E. Soper for 
correspondence. 
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Appendix A. Notations and Conventions in the Light-Cone Parton Model 

Four-vectors and scalar products are written in light-cone coordinates as 

1 



a , a, a ); a 



a- 



a 



+6+ + a h-+ a'bi = a+6 + a b'^ - a-b. 



(A.l) 



The Dirac 7-matrices are 



7 



7 



1 



V2 

1 

1 



(7° +7'); 7 



1 



7 



V2 
-a 
a' 



(7^^ - 7') 



(A.2) 



The Dirac field is decomposed into "good" and "bad" components by means 
of the projectors 

P+ = ^7"7+ P- = ^7+7"- (A.3) 



The good components ^_|_ and bad components ip- are then 






V^4/ 



^3 

V / 



(A.4) 



X 



+ 



Field operators are quantized in the light-cone gauge = 0, on the plane 
= (light- front quantization). Only the good components of the Dirac field are 
independent dynamical degrees of freedom. The bad components are given in terms 
of the good ones and the gauge degrees of freedom by the constraint equation 



1 



-i [{idj - gAj)'y^ + m] 7^+- 



(A.5) 



Color indices and factors, which are not relevant for our discussion, are omitted 
throughout. 

The basis of Dirac spinors appearing in the Fourier decomposition of 
Eq.(p?T|), is chosen to be 



/1\ 





Vo/ 



U- 






Vi/ 



(A.6) 



v± = u^, 

while fermion creation and annihilation operators are normalized according to 
'4(A;),4,(/c')} = {bsik),bl,ik')} = {2nf2k+5,s'S{k+ ~ k+')6(-''\k - k'). (A.7) 



20 



The basis of gluon polarization vectors is chosen as 



ei- = (o, 1, 0, ^) 
4 = (o, 0, 1, ^) . 
and gluon creation and annihilation operators are normalized as 



(A.8) 



(A.9) 
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TABLE CAPTION 



Table 1: The Feynman diagrams which contribute to Ba, Eq. (|3.18]) . The blob indi- 
cates insertion of the operator Bic, Eq.( |3.14|) , while the ordinary quark-gluon 
vertex is obtained from the insertion of the interaction Lagrangian Eq.( |3.19| ). 
± indicates the quark's helicity. The results listed correspond to the ampli- 
tudes {q\Cint\q9){q9\Bic\q) and {q\Bic\qg) {qg\Cint\q) , with Bic and Cint given 
by Eq.s ( |3.14| ) and ( |3.20| ), respectively. The states \q) are normalized according 
to \q) = b^s) with {s\s) = 1 in terms of creation operators normalized as in 

Eq.(EZD- 
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